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Abstract 

We provide a unified theory for the high force elasticity of biopolymers solely in terms of the 
persistence length, ^p, and the monomer spacing, a. When the force / > ~ ksT^^p/a^ the 
biopolymers behave as Freely Jointed Chains (FJCs) while in the range J^i ~ ksT/^p < f < Th 
the Worm-like Chain (WLC) is a better model. We show that can be estimated from the 
force extension curve (FEC) at the extension x ~ 1/2 (normalized by the contour length of the 
biopolymer). After validating the theory using simulations, we provide a quantitative analysis 
of the FECs for a diverse set of biopolymers (dsDNA, ssRNA, ssDNA, polysaccharides, and 
unstructured PEVK domain of titin) for x > 1/2. The success of a specific polymer model 
(FJC or WLC) to describe the FEC of a given biopolymer is naturally explained by the theory. 
Only by probing the response of biopolymers over a wide range of forces can the /-dependent 
elasticity be fully described. 

1 Introduction 

The elasticity of biopolymers, probed using single molecule methods, is often analyzed using the 
worm-like chain (WLC)^'^''^ model or the freely-jointed chain (FJC) model. For example, the force- 
extension curves (FECs) of double stranded (ds) DNA (at high monovalent salt concentration) are 
accurately reproduced using the WLC model with essentially the persistence length (^p) as the 
only adjustable parameter while the FJC model does not fit the data.= The WLC model is also 
successful in describing the FEC of ssRNA.^ and proteins.-'^ In contrast, the FEC of ssDNA is 
better explained using the FJC model than the WLC model. ^ Given that the values of ^j, of ssDNA 
and ssRNA are similar,— it is hard to understand the source of the discrepancy. In addition, for 
polysaccharide s iid2d3,i4 ^j^g FJC is used to extract the parameters from the FEC. 
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Here, we use general theoretical arguments to determine the polymer model that best describes 
the measured FEC. In particular, we derive a complete picture of the stretching physics of polymers 
at relative high forces, solely in terms of and the tensile screening length = ksT/f,^^ where 
ksT is the thermal energy and / is the stretching force. In the process, the FECs for the FJC and 
WLC models naturally appear as limiting cases within a unified theory for the stretching of polymer 
models, regardless of the nature of interactions between the monomers. Because the theory makes 
use of only general geometrical and scaling arguments that describe the response of a polymer under 
tension in terms of and ^t, we believe that our approach can be used to analyze the FECs of 
virtually any model polymer, and more importantly, many classes of biopolymers. 

We propose that the FEC for any polymer should exhibit a near universal behavior. At very 
small forces we expect the extension to scale as x ~ /.i£ii2ii^ As the force increases the FEC could 
exhibit a plateau region (poor solvent) or a universal behavior (good solvent) where the extension 
scales as X ^ j2/3i5 (^pj^cus regime). Here we show that this behavior cannot extend beyond the 
region x ^ 1/2. In the force range J-i ^ ksT/^p < f < Th ^ kBT^p/a^ with a being approximately 
the monomer size, the FEC should coincide with the prediction of the WLC model. For f > J-h 
the polymer is better described using the FJC model. Thus, as long as / > Ti, the FEC can be 
fit using the analytical results from WLC or FJC model regardless of the interactions that stabilize 
the ordered structures of the biopolymer at low forces {f < J-i). The model that provides the most 
accurate fit of the measured FEC depends on the force range of the FEC as well as the specific 
biopolymer, which is characterized by and a. In other words, the elasticity of the biopolymer (or 
any polymer for that matter) depends on the force.— 

2 Theory 

For a chain subject to a force / along the z-axis, the tensile energy is Hs = ^ fo-^iLi cos Of ^i, 
where 9f,i is the angle between bond i and the force vector. In the limit the tensile screening 
length = ksT/ f-i^^ is less than ^p, the excluded volume or other solvent- mediated intrachain 
interactions are irrelevant We make the physically reasonable assumption that, at high values 
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of /, the chain Hamiltonian has the global minimum when all the bond angles are zero, which 
corresponds to the rod state. For chains with special bond constraints such as the freely rotating 
chain, we shall assume that we can define virtual bonds and virtual bond angles that satisfy this 
assumption. 

2.1 The FJC regime is universal for all discrete chains at high force 

In the following we denote {X {e))ns = Jo X{e) exp sin OdO as the average with respect to 

TLs of X as an arbitrary function of some angle 9. For any two bonds i and j with bond angle 9ij, 
because of the triangle inequahty < % < df,i + dfj, the average of cos{9ij) with respect to only 
TLs satisfies: 

= (cos6'/.j}-Hs(cos6'/.j)-Hs - (sin6'/.j}-Hs(sin0/,i}-Hs 

and becomes independent of i and j. As the force becomes so large such that ^ a, Eq. ^ for 
any two bonds reduces to 

(cos0)«, >exp(^-^^^±il^). (2) 

On the other hand, the contributions due only to the intra-chain interactions lead to 

(cos6')-H,(s) = exp^-|-j, (3) 

where Tic is the chain Hamiltonian of any form and s ^ \i — j\a is the arc-length. The expression 
in Eq. ([3]), which defines the persistence length, is assumed to be true for many types of uncharged 
biopolymers and for polyelectrolytes at sufficiently high counter-ion concentrations. At high forces, 
the contribution from the tensile energy to the average in Eq. ([2]) becomes comparable to the 
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contribution from the chain interaction. Thus, the arc-length s at which the two contributions equal 
decreases as the force increases. Eventually, as the force approaches J-h such that s — a, we obtain 

exp (-^^^^) < {cos e)Hs {:Fh) - exp (^-^) , (4) 
which leads to an estimate of J-h , 

^. = c^, (5) 



with the constant c in the range from ^-t2^ to 4 (see below). For s > a, the average cosine is 
determined by the tensile energy Tis whose effects on the chain are (much) larger than intrachain 
interactions. At such high forces, the bonds align themselves along the force with small fluctuations. 
At forces above the stretching equation of a discrete semiflexible chain satisfying the above bond 
angle assumption becomes that of the FJC, i.e. 

x = (6) 
a 

Thus, should be viewed as a crossover force because above all biopolymers should behave as 
FJCs. Note that, in order for J-^ to be finite, a in equation ^ must not vanish, i.e., the polymer 
chain should be viewed as a series of connected discrete links. 

A similar form and notion of the crossover force were already pointed out elsewher e ^^'^^ for 
two specific models, WLC and freely-rotating chain (FRC). However, either the explicit form of the 
Hamiltonian (WLC) was required or simulations for FRC were needed to infer an empirical value of 
J-h- Here, general arguments based on the interplay between the force and intrinsic chain stiffness 
alone leads to an estimate of J-h (Eq.^). More importantly, J-h can be obtained without assuming 
any specific energy function for the polymer, and hence should be considered as a universal behavior 
for any discrete semiflexible chain. 

The critical value J-h defined in eqn. ([5]) could, in principle, be somewhat more complicated 
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especially if the bonds connecting the monomers are extensible. In this case, the stretching equation 
at high / for the FJC (Eq.®) has the formi^i^iSiiSS 




K' 



(7) 



where K is the stretch modulus that has the dimension of force. 

2.2 Force-dependent Kuhn segments and the universal effective FJC regime 

Consider the regime ^ < < (the lower bound is when / « Th)- Similar to the Pincus 
argument at low force (^t > Ik - the Kuhn length) in which the chain can be viewed as a series of 
blobs, the polymer breaks up into a series of segments whose length is on the order of the Kuhn 
length and behaves almost similarly to a FJC of Kuhn segments (see Fig.[T|). The relative extension 
in this force regime is, at most, a; ~ 1 — S^t/^K- However, we are not in the true FJC regime discussed 
above. The intrinsic intra-chain interactions, which tend to prevent the chain from bending, are 
relevant. The relative extension x must be larger than that of the FJC with bond length a, namely 
1 — ^t/a. Thus, we have the bound, 1 — ^<a;<l — In analogy with Eq. we can formally 



where A(/) is the segment length of the effective FJC, which should be regarded as an effective 
force-dependent Kuhn length. 

If the chain is viewed as a sequence of connected segments with effective length A(/) that re- 
sembles the FJC, the average cosine of the angle formed by any two of these A-segments is given by 
Eq. ([2]) where a is replaced by A: 



write 



X = 1 



A(/)' 



(8) 




(9) 
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The assumption that i^t <C A has to be justified a posteriori. The intrinsic tangent-tangent correlation 
of any two i and j efi^ective monomer-like A-segments is 



By equating the two equations ([9]) and (fTO|l (with \i — j\ = 1) we can relate A, / and using: 



It follows from Eg. pTj) that the force-dependent Kuhn segment Xk is determined not only by the 
intrinsic stiffness of the chain through but also by the tension through ■ The tensile screening 
length, which satisfies ft(oc /~^) <C A(/) (oc f^^^^) is similar to the deflection length proposed by 
Odijk— using the explicit quadratic form of the Hamiltonian. The FEC in this regime satisfies, 



which coincides with the results for the WLC when c = 4.— i^^i^^ Thus, by merely using the analogy to 
a semiflexible chain, which behaves as a FJC at high force and considering the effects of the intrinsic 
segment correlation and the effects of tension, we have obtained the characteristic stretching equation 
for any semiflexible chain. Such a regime exists in all semiflexible chains that satisfy the bond angle 
assumption. We shall refer to this as the effective FJC or WLC regime, generalizing the traditional 
name for any semiflexible chain. In practical applications the terms effective FJC and WLC can be 
used interchangeably. 




(10) 



(11) 




(12) 
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2.3 The unified stretching equation for both the FJC and WLC 

For the two regimes FJC and WLC discussed above, we assume that the generahzed equation for A 
takes the foUowing simple form, 

A(/,a) « (A^+a")^, (13) 

where a > 1 is a model-dependent parameter. This equation reflects the observation that when 
/ < J-in A(/) is almost Xk a) and the chain behaves as an effective FJC; and when / > Th, 
A(/) tends to a {':^ Xk) and the chain behaves as a true FJC. Thus, we have the equation for the 
two strong stretching regimes: 

x^l (14) 

(A^ + a")° 

which reduces exactly to the form for the discrete WLG^S when a = 2. 

Because of the upper limit of A, which is Ik = 2^p, the lower limit on the force for observing the 
effective FJC regime (assuming c < 4) is: 

X^lK^^t^^p^J'i^ (15) 

Here, can be called the crossover force to the effective FJC regime from a low force regime (see 
Fig.[l]). 

The physics of the stretching of all biopolymer molecules using discrete semiflexible chains is 
summarized in Fig[TJ At very low forces, the extension increases linearly with the force. For larger 
forces, but not too large, such that the tensile screening length is still significantly larger than the 
Kuhn length, one reaches the Pincus regime of blobs where the extension scales either linearly with 
the force (ideal chains) or as f^^^ (self-avoiding chains).— As the force increases until it is larger 
than Ti, becomes smaller than the Kuhn length, the chain behaves as an effective FJC with the 
law 1 — a; oc /~^^^. Eventually, when S,t becomes smaller than a^/^p (/ > Th), the chain behaves as 
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a true FJC under high tension with the FEC law 1 — x cc f 



3 Simulations confirm the theoretical predictions 

To test our theory, we performed simulations using a number of models. 

"Polynomial" models ("Generalized" WLC). The Hamiltonian of the chain with N bonds is 



1=1 



where 9, the angle between two successive bonds, varies from to tt. When 7 = 2 and 8 = 1, the 
Hamiltonian approximately reduces to the WLC energy function and when 7 is large, we recover 



the local approximation Thick Chain model.— Fig. 2(a) shows the dependence of Tib on 6/Q for 
several values of 7. 

The FBNE (finite bendable non-linear elastic) chain. To model the finite bendability of a chain, 
we use the FBNE potential 

00, o>e, 

which resembles the FENE potential— 1^" for chain bonds. 

To compare all simulation data, we chose parameters in the models such that the persistence 
length is fixed at £_p/a ~ 10. The ratio ThjJ-i ~ c{^p/a)'^ = 100c is large enough to clearly 



observe the effective FJC regime. From equations (|12|) and ([6]), it can be seen that plotting 1 — x 
against / x a/ksT in double logarithmic scale would give us the characteristic slopes of —1/2 and 



— 1 for the effective FJC and FJC regimes, respectively. Indeed, Figs. 2(c) and 2(d) shows that 
there are these characteristic slopes for the "polynomial" model at all 7 considered. The same 
result is obtained for the FBNE and self-avoiding TC models^'^'^^ (data not shown). Moreover, 
the fits using the generalized equation (|14p with two free parameters c and a are excellent and 
yield, in accord with theoretical arguments, c in the range to 4, and a in the range 2 to 7 
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and positively correlated with 7. Moreover, it is easily seen that the effective FJC regime sets 
in as soon as f ~ 0.1-^^^ ~ a_kBT_ ^ -p ^\^^ crossover from effective FJC to FJC occurs at 
f ~ 40-^^ ^ _^K^BJ_ ^ jr Thus, the simulations confirm the theoretical arguments. 

4 Analysis of experiments 

4.1 Significance of the regimes and crossover forces 

Our analysis provides a theoretical basis for choosing the appropriate polymer model to analyze 
the FECs of biopolymers. In order to choose a specific model it is necessary to estimate and 
J-fi, which require the values of a and (see eqs. (jlSp and ([5])) that are intrinsic properties of the 
biopolymer. The value of the monomer spacing a can be obtained by analyzing the FEC at high 
force (FJC regime). The effects of the self-avoidance of the chain can be assessed most accurately in 
the linear and Pincus regimes (FigU]) . Chain extensibility can make it difficult to estimate a making 
it necessary to include models that treat chain extensibility explicitly.— iSli^ The persistence length 
is easily extracted from the FECs in the effective FJC regime. It is the most accurate quantity that 
can be inferred from experiments, since it is the effective FJC regime that is commonly detectable in 
those experiments.— i^'^'^i^iii^'^i^^ In practice, we use the force at a; « 0.5 (referred to as 1/2-rule) 
to estimate Ti ~ ksT/^p^ and hence (see Ec^. (IT^ with c w 4). We illustrate the utility of our 
theory that depends on the interplay of the two "crossover" forces Ti and Th in determining the 
elasticity of biopolymers by analyzing measured FECs. 

dsDNA: Using the well-known values « 50nmiii^iM and a = 0.34 nm at T 300K, we 
obtain Ti « O.lpN and Th ~ 6400pN. Values of / that are comparable to J-i are accessible only 
using magnetic beads or optical tweezers, ^i^i^ whereas / ^ J-h can be realized in AFM experiments. 
Because dsDNA undergoes a phase transition to an overstretching form at / « 65 pN,— i^^i^i^^ the 
FJC limit cannot be observed in dsDNA. It should be stressed that even if a is greatly increased, up 
to the conceivably maximum value of 3.4 nm or one double helix turn, Th still exceeds / ~ 65pN. 
These estimates show that the FEC of dsDNA (excluding the overstretching region) should be in 
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only the effective FJC regime, wliicli explains the success of the WLC in quantitatively describing 
the FEC of dsDNA. 

ssRNA: From the typical values of S^p « Inm— at sufficiently high counter- ion concentrations 
and a « 0.55 nmi^ for both ssRNA and ssDNA we obtain J^; « 4 pN and ~ 50 pN. Thus, 
both J-'i and J-h are small enough that they can be accessed in typical LOT experiments, making it 
practical to infer the /-dependent changes in the elasticity. The values of J-i and J-h readily explain 
why the measured FEC for poly-U using optical tweezers^- in the force range from about 0.5 pN 
to 50 pN, which is about Th, is best described using the WLC model. We predict that the FJC 
behavior should emerge for forces larger than 50 pN for poly-U. The same arguments also hold good 
for ssDNA data at high ionic concentrations in the following paper by Saleh et al. 

ssDNA: In a recent experiment, Saleh et al. measured FEC of ssDNA at varying counter-ion 
concentrations in the force range 0.1 ^ f ^ 50pN.— The data suggests that for / < /c, a critical 
force that depends on the salt concentration, the FEC is in the Pincus regime (Fig. [1]) in which the 
non-linear response to force is determined largely by the excluded-volume interactions. We adopt 
the 1/2-rule to calculate ^p{Xd) as a function of the Debye length Ad for ssDNA and poly-U. The 



estimated ^p(Ad) for these two polyelectrolytes are identical (Fig. 3(a)). Moreover, we can fit the 
curves with ^p{\d) ~ ^p+^^D, where is the "bare" persistence length and A is a proportionality 
constant. The values of A and for poly-U and ssDNA are nearly identical (^^'^"'''"^ — 0.67nm and 
^o,ssDNA _ 0.63nm). The linear dependence of £_p on Ad is in agreement with previous works^Sdl as 
well as Saleh's findings.-Sfl Although the procedure used by Saleh et al. to extract £,p{\d) is different 
from the WLC-based analysis presented here, our estimated ^^^''^'DNA -^^ excellent agreement with 
their estimate of 0.6nm.— 

Comparison of our predictions with the FEC data for ssDNA reported by Rief et al. ''^ in the force 
range from lOpN to 200pN and in Tris buffer further validates the present analysis. In a different 
paper ,~ Dessinges et al. combined the data with their data on ssDNA in the force range from 0.05pN 
to 20pN and in a different buffer (see Fig. 4 in ref.— ). Using extensive data on ssDNA, they showed 
that at very low force, self-avoidance and electrostatic effects play dominant roles, and hence the 
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elasticity of ssDNA is qualitatively different from the high force regimes. Indeed, by taking into 
account the electrostatic interactions in the so-called modified (extensible) FJC model, Dessinges et 
al. explained the low force part of the FEC. They also showed that both the WLC and the "bare" 
modified FJC, when applied independently, failed to fit the entire high force part {x > 1/2) of 
the FEC. Here we show that the FEC above f = Ti (data by Rief et al.— ) can be quantitatively 
analyzed using our theory. The equation for the high stretch regimes (Eq. (|14p ). which has both 
the WLC and FJC high force behaviors, provides a very good fit to the data. We notice that the 
modified FJC model by Dessinges et al. with electrostatic interactions actually failed to fit the high 
force {x > 1/2) portion of the FEC* It could be argued that the good fit to the lower portion of 
the FEC was mainly because the electrostatic interactions, which play dominant roles, were well 
taken into account. The fact that whether the FJC or WLC is used for the low force regime should 
not considerably affect the outcome of the fit. But for the high force regime, it is essential to do a 
force range analysis prior to the actual fit, as it is done here. Our extracted values are a « 0.5 nm, 
£,p ~ 0.8 nm and a w 7 leading to J^; « 5 pN and ~ 55 pN (see Fig |3(b)[ ). 

Protein (PEVK domain): PEVK (Pro-Glu-Val-Lys) domain is a proline-rich domain of 
the muscle protein titin and it has no definite folded structure. Using the accepted values « 
0.5nm— iii^^i^ and a = 0.38 nm, « 8.4 pN, which is well in the covered range of almost all single 
molecule apparatuses, and w 50 pN. We expect that the WLC is more suitable for the PEVK 
domain than the FJC, and indeed FECs in the force range from 10 to 200pN were analyzed using 
WLC.-^'i'^'^a However, we can see the need for a discrete model rather than the continuum WLC 
to obtain a better fit. 

Polysaccharides: Linear polysaccharides that occur naturally in a variety of cellular structures 
are subject to tensile stress. At high forces the sugar rings undergo a transition from chair to boat 
conformation, which clearly changes their elastic properties.— Unlike other systems (dsDNA, for 
example) the change in elasticity is associated with enthalpic changes as a result of the conformational 
transition. Surprisingly, it is found that over a wide range of forces the FEC is best explained by 
the FJC model, which can be fully explained using our theory. For cellulose with the value 
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Ik K, a K, 0.54nm,— liSiii^iii we obtain J^i « J^h ~ 15.3pN. This means that excluding the overstretch 
range, the entire range of FEC for cellulose reported in references by Marszalek et al.— li^ii^iii^ 
obtained using AFM with force range of 10 to 1000 pN is best described using the FJC. The same 
observation can be made for other polysaccharides, such as dextran {Ik ~ a « 0.44 nm in the chair 
conformation and Ik ~ a ~ 0.57 nm in the boat conformation) and amylose [Ik ~ a ~ 0.45 nm 
(chair) and Ik ~ a~ 0.54 nm (boat)).!- Thus, the intrinsic properties of polysaccharides show that 
this class of molecules behave as FJCs. 

In all the above examples, it is necessary to assess the range of force in any particular FEC 
before using the polymer model that can best describe the data. In most cases, it is necessary to 
consider both the semi-flexibility and discreteness of the biopolymer. Models such as the discrete 
WLC modeP^ and the TC model— are better for most fitting purposes. In analyzing most of 
the currently available data we find that these two models are adequate. However, as shown by 
Dessinges et al.— the variations in the elasticity of ssDNA over a wide range of forces will require a 
combination of models. Our theory accounts for /-dependent variations in the elasticity in terms of 
naturally occurring length scales as long as a; > 1/2. Below x = 1/2, additional forces may stabilize 
specific structures. 

4.2 Overstretching in a FEC 

For each biopolymer, there is a threshold force at which the transition to bond overstretching occurs. 
Inextensible models can be used to fit only the experimental FECs below the threshold force. From 
the approximate equation in the high stretch regime it can be shown that the slope of the 
FEC above J-i in the log-linear scale is a strictly monotonically increasing quantity for a generic 
semiflexible chain. On the other hand, in the region where overstretching effect is significant, the 
extension can be approximated as^ x ~ 1 — jtt^ + "^i where K, m and C are positive constants. 
It can be shown that the rate of change of the slope in the region a; > 1 (overstretched) is always 
negative. Then for a FEC with overstretching, the threshold force is the one (> J^i) at which the 
slope of the curve In /(a;) vs. x stops increasing. A slope that stays constant or decreases would 
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signify the failure of the inextensible model. 

The simple analysis, that accounts for extensibility through the stretch modulus allows us to 
gauge the threshold force from the FEC. If the data are precise, we can take the numerical derivative 
of In / (a;) in the region f > J-i and detect the force at which it reaches a maximum value or a plateau. 
Otherwise, the practical solution would be to plot f{x) vs x in the log- linear plane then detect the 
change in the slope above Ti by inspection. Although not rigorous, it turns out to be an effective 
way as demonstrated in Figured] For the poly-U data^ at various salt concentrations and at forces 
/ < 50pN, , the practical technique reveals no decrease in the slope in the x > 1/2 region (data not 
shown). Therefore, we expect no overstretching in poly-U or ssDNA. 

5 Conclusions 

Our theory shows that the WLC and FJC models at large forces are the two limiting universal 
behaviors for virtually all polymers and especially biopolymcrs. The crossover from WLC to FJC is 
determined by the persistence length and the monomer spacing a. The validity of the theory is 
strongly demonstrated through the excellent agreement with simulation results and the successful 
analysis of the force extension curves of a diverse set of biopolymers. The two crossover forces, Th 
and Ti^ which demarcate the boundary between the applicability of the WLC and FJC, provide a 
physical picture of the /-dependent elasticity of biopolymers. The derived 1/2-rule is a convenient 
way to obtain a good estimate of the persistence length of the biopolymer from the measured FEC, 
which in turn yields the model that is most appropriate for analyzing the FEC. The theory also 
shows that complete understanding of the elasticity of biopolymers requires measurement of FECs 
for a wide range of forces as obtained for ssDNA.— 
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Figure captions 



Figure [T] Summary of the various stretching regimes in polymers or biopolymers. In the smallest 
force regime, the extension increases linearly with the force. In the case of biopolymers (RNA and 
protein) the structure in the collapsed state is ordered. Polymers in poor solvents are collapsed. 
The size of the polymer or a globular protein is about N'^a, where N is the number of monomers, 
a is the monomer spacing and v is the scaling exponent, which is w 1/3 in poor solvent, 1/2 in 
6'-solvent or for ideal chains, and 3/5 in good solvent. In good or 6'-solvent, as the force increases 
the chain extends into a series of Pincus blobs with typical size — ksT/ f. In this regime, the 
extension scales as a; oc ji/'^-i, i5,i7 \Yhen the force is comparable to J-i = ksT/^p {^t ~ (.p), each 
Pincus "blob" contains just a single segment on the order of a persistence length. In this force range, 
the Pincus scaling law breaks down and the elasticity if described by the WLC model. The chain 
can be viewed as an effective freely jointed chain of segments of length Xk ^ \/it£,p- Intrachain 
interactions may still be relevant within the length scale Xk but are negligible on larger length 
scales. The extension goes as 1 — x oc f^^^^- Finally, when the force reaches J-h ~ ksTS^p/ a? ^ all 
chain bonds are strongly aligned with the force and the intrinsic chain interactions are irrelevant. 
Thus, the chain behaves as a freely jointed chain (segment length a) with the extension scaling as 
1 — X (X- f^^ ■ In some cases of biopolymers discussed in the main text, Tu can be so high such that 
bond overstretching can occur before the force reaches Th and the theoretical FJC regime may not 
be actually observed. 

Figure [2j Simulation of model Hamiltonians and fits to the simulation data, (a) Plot of 
Hb/Ki, as a function of 9/Q at various values of 7 (Eq. (fT6|) ). As 7 increases from 2, which is 
approximately the value for the WLC model, the potential mimics the excluded volume interaction 
between consecutive monomers, (b) Plot of 'HFBNE{0)/kBT in Eq. fTTl) at two sets of values of Kt 
and Q with £,p/a = 10. (c) Log-log plot of 1 — a; as a function of fa/ksT obtained from Monte 
Carlo simulations for the polynomial model with ^p/a = 10 and — 400 and various values of 7. 
The scaling laws 1 ~ x /"^^^ and 1 — x ^ /^^7 as well as the two crossover forces Ti and are 
clearly visible, (d) The fits of Eq. (fT4| to the simulation data for 7 = 2 and 20. 

Figure [3l Analysis of the experimental results, (a) Estimates oi£,p{Xo) for poly-U (circles) and 
ssDNA (filled symbols) using the 1/2-rule. Linear fits to the data (solid and dashed lines) yielded 
the bare persistence lengths ^o>P°iy-u « 0.67 nm and ^O'S^^dna ^ o.63. (b) The fit of equation ^ 
(solid line) to the EEC of a ssDNA (squares) by Rief et al.--The extracted values for the parameters 
are: « 0.8nm, a ~ 0.5 nm and a 7 leading to ^/ « 5 pN and J-'h ~ 55 pN. The extracted values 
of S.P and a are consistent with those reported in ref.JS and,"^^ respectively. The value of a ~ 7 
implies that 7 (see Eq. 1161) is large. As a result, the bending energy increases more sharply than 
predicted by the WLC model (Fig. 2(a) ) as 6 increases beyond a critical value. 

Figure [4l An illustration of overstretching in cellulose. A threshold is identified as a narrow 
region in the circle, where it shows a decrease in the slope of the FEC when plotted in log-linear 
scale. 
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